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ABSTRACT 

We model the time-dependent evolution of metal-enriched intergalactic and circum- 
galactic gas exposed to the fluctuating radiation field from an active galactic nucleus 
(AGN). We consider diffuse gas densities (nn = 10^'"' — lO^-^cm^-^) exposed to the 
extra-galactic background (EGB) and initially in thermal equilibrium (T ^ 10^ — 10*'^ 
K) . Once the proximate AGN field turns on, additional photo- ionisation rapidly ionises 
the H I and metals. The enhanced AGN radiation field turns off after a typical AGN 
lifetime (tagn = 1 — 20 Myr) and the field returns to the EGB intensity, but the metals 
remain out of ionisation equilibrium for timescales that can significantly exceed tagn • 
We define this phase as the AGN proximity zone "fossil" phase and show that high 
ionisation stages (e.g. Ovi, Neviii, Mgx) are in general enhanced, while the abun- 
dances of low ions (e.g. Civ, Oiv, Mgii) are reduced. In contrast, Hi re-equilibrates 
rapidly (<C tagn) owing to its low neutral fraction at diffuse densities. We demon- 
strate that metal column densities of intervening gas observed in absorption in quasar 
sight lines are significantly affected by delayed recombination for a wide range of den- 
sities, metallicities, AGN strengths, AGN lifetimes, and AGN duty cycles. We mo del 
the exceptionally strong z = 0.9 Neviii absorbers observed bv iTripp et al.l (|201lll as 
arising in a possible fossil zone or near a recently turned-on AGN and we demonstrate 
that at low redshift even moderate-strength AGN that are off for 90% of the time 
could significantly enhance the high-ion metal columns in the circum-galactic media 
of galaxies observed without active AGN. Fossil proximity zones may be particularly 
important during the quasar era, z ~ 2 — 5. Indeed, we demonstrate that at these 
redshifts a large fraction of the metal-enriched intergalactic medium may consist of 
out-of-equilibrium fossil zones. AGN proximity zone fossils allow a whole new class 
of non-equilibrium solutions that may be applicable to a large fraction of observed 
metal absorbers and which could potentially change the inferred physical conditions 
and masses of diffuse gases. 

Key words: atomic processes; plasmas; galaxies: formation; intergalactic medium; 
quasars: absorption lines; Seyfert; cosmology: theory; 



1 INTRODUCTION 

Active galactic nuclei (AGN) provide the majority of 
photons capable of ionising intergalact ic hydrogen over 
most of the age of th e Universe (e.g. IShull et al] 1 19991 : 
iHaardt fc Madaull2012l ). Some of the most commonly ob- 
served metal absorption lines in the intergalactic medium 
(IGM), including C iv and O vi, have ionisation potentials at 
extreme UV (EUV) energies where their main photo-ionising 
sources are AGN. The mean free paths of ionising pho- 
tons with energies greater than 1 Rydberg are long enough 
to warrant the assumption that the ionising extra-galactic 
backgrou nd (EGB) is uniform in its spatial distribution (e.g. 
IZuoll992l ). The main criterion for this assumption is that the 



mean free paths of photons are lon ger than the clustering 
lengt h of photo-ionising sources (e.g. iFaucher-Giguere et "all 
l2009h . 

This assumption obviously breaks down near sources 
of ionising radiation. Particularly notable are the proxim- 
ity zones of QSOs, where the locally generated ionising 
radiation at EUV energies may exceed the EGB by or- 
ders of magnitude. The simple nature of hydrogen, hav- 
ing only two ionisation states, allows an elegant relation- 
ship where, in ionisation equilibrium, the neutral fraction 
scales in inverse linear proportion to the ionisation field 
strength if hydrogen is photo-ionised. A dramatic example 
of this are the proximity zones of high-z quasars, where 
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the Lya forest absorption is significantly reduced by the 
enha nced AGN radiation field as the QSO is approached 
(e.g. ICarswell et al.l 1 19871: iBaitlik et all Il988l : IScott et alj 
I2OO0I : iBolton fc Haehneltll2007l ). QSO pairs can probe the 
transverse proximity effect in which the background QSO 
line of sight intersects the foreground QSO proximity zone, 
where the radiation field is increased by orders of magnitude 



phase, which is often calcu l ated to be 



10" 



10** 



compared to the typical EGB (e.g . iJakobsen et alJ 20031: 
Schirber et al] |2004| : iHennawi et al]|2006l : iGoncalves erall 



200^ ). 



Metal absorption also reflects the effects of increased 
ionisation in quasar proximity zones. Proximate absorbers 
are often defined as being within 5000 kms~^ of the QSO 
velocity along the line of sight, and show some distinct sig- 
natures s uch as O VI absorption with little or no H i absorp- 
tion (e.g. lTripp et al.ll2008l ). IGoncalves et al.l (|2008l ') observe 
unique metal-line systems at z ~ 2.5 in transverse proximate 
zones with an ionisation enhancement of 10-200 x, where 
O VI is enhanced and C iv is weakened. These signatures 
suggest that enhanced photo-ionisation affects metal lines, 
but in a different way than hydrogen. The multiple ionisa- 
tion levels of heavy elements makes the response to increased 
ionisation more complex: metal ion fractions of order unity 
are possible at diffuse IGM densities (i.e. much higher than 
H I fractions) ; plus these fractions are not a monotonic func- 
tion of density or temperature like that of H I, but instead 
peak at a preferred density or temperature for an interme- 
diate ionisation state. In the case of an increased ionisation 
field, the peak ionisation fraction for a metal species will 
move in linear proportion with the photo-ionisation rate to 
higher density if the species is photo-ionised and in equilib- 
rium. 

Another difference between hydrogen and metal species 
is the timescale that it takes each to achieve ionisation equi- 
librium when a proximate AGN turns on and off. The ioni- 
sation timescales, relevant for the AGN turn on, are longer 
for metal species with higher ionisation potentials, because 
fewer ionising photons are produced by the AGN at these 
higher energies. However, the more important difference is 
the longer timescales metals take to reach ionisation equilib- 
rium after the AGN turns off. Although most metal species 
have shorter recombination timescales. Tree, than Hll, the 
relevant timescale to achieve equilibrium is Tree times the 
equilibrium ionisation fraction of the recombined species. 
For example, at typical warm (T ~ 10" K) IGM densities 
of the Lya forest (nn ~ 10~*cm^''), Tree, hi: ~ 10® yr, but 
the Hi ionisation fraction, /hi, is ~ lO"**; hence the equili- 
bration timescale is only ~ 10'' yr, because hydrogen only 
has to recombine to a very low neutral fraction to achieve 
equilibrium. Metal-line species like C iv and O vi have re- 
combination timescales that are > 10^ yr at similar warm 
IGM densities, but typical ionisation fractions much closer 
to unity. Therefore, timescales to recombine to equilibrium 
are also > 10^ yr - much longer than for Hi. Furthermore, 
the existence of multiple ionisation states for metals means 
that there are a series of ions to recombine through, and 
hence multiple recombination timescales, further extending 
the total equilibration timescale. 

AGN proximity zone fossils are defined as proxim- 
ity zones in which the AGN has turned off, but ioni- 
sation equilibrium has not yet been achieved. The fos- 
sil lifetime can last much longer than the AGN-on 



yr (e.g. iHaiman fc Hull I2OOII: iMartini fc Weinberg! 12001 



Jakobsen et all l2003l: ISchirber et all 12004 iHopkins et all 
20061 : iGoncai ves et al.l l2008h . In th s work, we use the 
non-eq uilibrium code introduced in lOppenheimer fc Schavd 
(|2013t ) to follow the ionisation and temperature evolution of 
metal-enriched gas subjected to a variable ionisation field 
owing to an AGN. We consider several situations where the 
fossil stage lasts as long as, and sometimes much longer 
than the AGN phase. We predict the absorption line sig- 
natures of these fossil zones, and compare them to some 
recent observational results. We focus on diffuse gas with 



hydrogen numbers densities nu — 10 



10" 



(i.e. 



p/{p) ~ (50-5000) X (l-l-z)"^), initially at the equilibrium 
temperature where cooling balances photo-heating, corre- 
sponding to warm IGM temperatures, T ~ lO** — lO'*'^ K. 

The first case we will consider may be applicable to 
the case of the recently discov ered strong Nevill absorber 
detected bv lTripp et all l|201lh in the Cosmic Origins Spec- 
trograph (COS) spectrum of PG1206-I-459 at z ~ 0.927. 
This system is one of the most unique and most puzzling 
metal absorbers ever observed, containing both high ions like 
Nevill and low ions like Mgll. We explore non-equilibrium 
photo-ionised models with the AGN on and in the fossil 
phase. We bolster our argument with evidence for a post- 
starburst, > L* galaxy, possibly within 100 proper kpc, 
which shows signatures of being an AGN. We use this as 
a fiducial example to consider variations in AGN strength, 
lifetime, and duty cycles to show how metal lines can re- 
spond to variable ionisation conditions. 

We next consider AGN proximity zone fossils at the 
height of QSO activity at 2 = 2.5. The observations of metal- 
line absorption in transverse QSO proximity zones along the 



line of sight (e.g. Worseck fc Wisotzkil 120061 : IWorseck et al.l 



I2OO7I : iGoncai ves et al.ll2008l ) provide our initial motivation, 
and we extend the exploration to the case where the prox- 
imity zones turn into fossil zones after the foreground QSOs 
turn off. We argue that unique absorption signatures should 
occur where Hi re-equilibrates, but Ovi and Civ are still 
out of equilibrium. We suggest that during the AGN era a 
large fraction of metal-line systems could reside in proximity 
zone fossils. 

The last example we consider concerns massive, spiral 
galaxies at z = 0.25 that may recently have harboured active 
AGN, such as Seyferts, but are currently in the less active or 
off phase. We demonstrate that increases of at least a dex 
are possible for circumgalactic (CGM) O vi time-averaged 
over fossil phases for normal Seyfert luminosities turning 
on 1 Myr out of every 10 Myr. Our examples show just 
how sensitive metal lines can be to moderate AGN, even if 
they are only on for a small fraction of time. As nearly all 
galaxies are thought to contain supermassive black holes in 
their centres, fossil proximity zones may be very important 
for the interpretation of QSO absorption line data. 

This paper is organised as follows. We introduce our 
non-equilibrium code and the methods used in this paper 
in ij2] We apply these methods to the AGN proximity zones 
and proximity zone fossils for the examples mentioned above 
in ij3] and we summarise in Sjl] Throughout we use proper 
distances unless noted otherwise. 
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2 METHOD 

We use the n oii-cquilibrium solver introduced in 
lOppenheimer fc~Sc havG ( 2013) to follow the time-dependent 
ionisation states of 11 elements (H, He, C, N, O, Ne, Mg, 
Si, S, Ca, and Fe). The solver also calculates ion- by- ion 
cooling and photo-heating to follow the time-dependent 
non-equilibrium net cooling, defined as cooling minus 
photo-heating. The code is intended to be integrated in 
hydrodynamic codes, and includes a treatment to sub-cycle 
the integration of ionisation and cooling given the input 
Courant timestep. 

In this work, we add the functionality to our method 
to allow the addition of photo-ionisation and photo-heating 
from a time variable source, which in this case has the 
spectral energy distribution (SED) of an AGN. The AGN 
photo-ionisation rate (s~^) for an ionisation state i of atomic 
species x from a radiation source with specific intensity 
(erg s~^ cm~^ Hz~^) is 



f°° fv 



(1) 



and the AGN photo- heating rate (erg s ^) is 
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(2) 



where v is frequency, i'o,xi is the ionisation frequency, ^^^(i^) 
is photo-ionisation cross-section, and h is the Planck con- 
stant. Since we track 133 ionisation states of 11 elements, 
we integrate the above equations for the 122 ionisation states 
with electrons. 

For the EGB, we replace with the ionising back- 
ground SED averaged over the whole sky, iirj^, in 



Equations [T] and [21 where J, 
-1 



has units of erg s cm 



Hz 



sr As explained in lOppenheimer fc Schavd (|2013h . 
we calculate Pegs and eBOB at 50 different redshifts 
between z = ~ 9 using the iHaardt fc Madaul (|200ll . 
hereafter HMOl) EGB and the s ame a^i (v) used i n 
CLOUDY ver. 10.00 last described in lFerland et all (|l998f ). 
We use the HMOl EGB as our fiducial background since 
this results in good agreem ent with the obse r ved column 
density distribution of H I (|Dave et al.l l2010l : Altav et al.l 
2011 ) and metal line streng t hs an d ratios JSchaye et al.l 



20031: lOppenheimer fc Davd l2006l: lAguirre et al 



20091 : lOppenheimer et al.l 



2008 



2012 



Oppcnhcimer & Davej 
Rahmati ct al. 20ii). 

The time dependent evolution of the particle number 
density, n, for each ionisation state, Xi, is given by 



duxi 



dt 

+ Pa;i_i,AGN) — U^i {{Oixi + Px^ris + Pa;-, EGB -|- Pa;;, AGN ) 

(3) 

where rie is the free electron density (cm~^), a^,. is the to- 
tal recombination rate coefficient (radiative plus di-electric, 
cm'^ s~^), and is the collisional ionisation rate coef- 
ficient (cm^ s~^). The AGN subscript indicates that the 
coefRcients have been calculated using our AGN template 
spectrum explained below. Gharge transfer and Auger ion- 



Figure 1. The lHaardt fc Madaul 1 I2OOII) model for the spectrum 
of the EGB closest to the three redshift regimes (z = 2.5, 0.9 
& 0.25) we consider, along with our AGN template spectrum, 
normalised to Jll = 10~-^^ ''ergs~^ cm~^ Hz~^ sr~^. The EUV 
power law of the AGN i s between 1300 A and 354 eV 

based on observations bv lTelfer et al ] |2003). The X-ray portion 
of the spectrum is assumed to follow the CLOUDY AGN template 
(j'^" '^ until 100 keV and then i/^^'^'' at higher energies). The 
four coloured vertical lines at the bottom indicate the ionisation 
potentials needed to reach the corresponding ion listed. 



isa tion are also included in our calculations as explained 
in lOppenheimer fc Schavd (|2013l ). but were omitted from 
Equation |3] for simplicity. 

The net cooling rate per unit volume (Anet, erg cm""^ 
s~^), is the difference between cooling and photo-heating 
from the EGB and the AGN, 



J^nct,Xi{T, z,nx-,ne) = A'^.{T)nx-ne — exi,EGB{z)nx- 
— ^Xi,AGN{z)nx.. 



(4) 



The units of the cooling efficiency per ion, A^. are erg cm"^ 
s~^, where we divided the cooling rate per volume (erg cm~'^ 
s~^) for ion Xi by nx-Ue to achieve density independence of 
collisional cooling, which allows us to tabulate these rates 
as a function of temperature only. 

Figure [T] compares our AGN template spectrum, 
to the HMOl EGB, J„{z) at the three different redshifts 
{z = 2.48, 0.87, and 0.27) nearest the three redshifts we 
consider in this paper {z = 2.5, 0.9, and 0.25). The HMOl 
EGB strength at the Lyman limit at these three redshifts 
is from high to low-z log[JLL (ergs~^ cm~^ Hz"'^ sr~'^)] = 
-21.30,-21.64, and -22.16. 

In cyan is the AGN SED normalised to /ll ~ 
■|^g-20.9gj.gg-i |,jjj-2 jj^-i^ Since this is a point source, 

we can convert its flux into Jll equivalent units 
by dividing by An steradians to get JLL.cqu ~ 
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■[^Q-22.0gj.g g-i jj2~^ sr^^, which has the same units as 

the EGB. Hence, JLL.equ is the intensity the EGB would 
have if it were equaUy bright as the AGN at the Lyman 
hmit. The spectrum has an EUV power law with slope 
^- 1.57 |-)g^ygg]-^ 1300 A and 354 eV based on observations 
bv lTelfer et"al] (|2002l ). At higher energies, the spectrum fol- 
lows until 100 keV and then at higher energies, 
which follows the template CLOUDY AGN spectrum. The 
most important slope for our consideration is the EUV (i.e. 
10-354 eV) slope, because these energies correspond to the 
ionisation potentials of the metal species we consider and 
which are marked in Figure [T] However, the X-ray portion 
of the spectrum is also important for two reasons: i) photo- 
heating rates of species like O vil and O vill depend on the 
shape of the X-ray spectrum, and ii) we convert to observed 
luminosities of AGN based on their 0.5-2.0 keV luminosities 

{Lo.5-2.o)- 

To convert these units to AGN luminosities, the above 
flux (/ll = 10~^'''^erg s~^ cm~^ Hz~^) corresponds to 
io. 5-2.0 = 10^^'^ergs^^ AGN at 100 proper kpc, which is 
the luminosity one may expect from a local Seyfert. The 
Seyfert would outshine the z = 0.25 HMOl EGB at 912 
A by 45% at 100 kpc, but because of the harder AGN spec- 
trum at the Ovi photo-ionisation edge, 114 eV, the AGN 
would outshine the HMOl spectrum by a factor of 23 x, as- 
suming no absorption of the AGN spectrum. In general, for 
our template AGN spectrum the conversion from AGN soft 
X-ray luminosity to equivalent EGB intensity units is 

T _ /ll _ ,p,-22 f TAGN \ ^ f Lq. 5-2.0 \ 

-^LL.e<5" - 4^ - y^^^ J 1^1042.iergs-i J 

-1 -2 -rj -1 -1 

ergs cm Hz sr , 

(5) 

where tagn is the proper distance from the AGN. 



3 APPLICATIONS OF AGN PROXIMITY 
ZONE FOSSILS 

We consider three situations spanning a wide range of red- 
shifts, from 2 — 2.5 to 0.25, that are applicable to recent 
observations. In all cases we assume the gas is being irradi- 
ated by the HMOl EGB for the given redshift, and that it 
initially has the equilibrium temperature for which photo- 
heating balances cooling. This temperature corresponds to 
T ~ ]^q3.85-4.68 assuming solar abundances for the densi- 
ties we consider, nn = 10~^" — 10~^"cm~''. The main rea- 
soning behind our parameter choices is that simulations indi- 
cate that a significant fraction of diffuse metals reside along 
a locus in density-temperature phase s pace where photo- 
heating balances cooling ( e.g. lOppenheimer fc Pavel [200^ : 
|w lersma et al.ll2010l . 120111 ). and that many commonly ob- 
served metal-line spec ies arise from gas at the se densities 
and temperatures (e.g. lOppenheimer et al.ll2012 ll. Although 
solar abundances are used, we will argue that our theoretical 
metal columns scale nearly linearly with metallicity. 

In i]3.1l we discuss our first application: an AGN turning 
on at 2 = 0.9, which is motivated by observations of the very 
strong Ne vill system observed bv I Tripp et al.l (|201ll ') toward 
PG1206-I-459, which appears to reside near a galaxy showing 
signs of AGN activity. We pay particular attention to this 



example, because of the large number of metal-line species 
observed in this system and also because Ne Vlll is one of the 
species most affected by AGN ionisation. We vary several pa- 
rameters including metallicity, density, AGN strength, and 
AGN lifetime. In ^3.21 we consider the effects on C IV, N v, 
and O VI columns at z — 2.5, both when a quasar is on and 
after it turns off. We estimate the fraction of the volume 
of the IGM altered by fossil zones, and argue that a large 
fraction of observed metal lines may be affected. Lastly, in 
i|3.3l we explore short duty cycles of 2 = 0.25 AGN that 
significantly enhance Ovi columns, even while the AGN is 
off. This case could be applicable to normal looking spirals 
that were recently Seyferts, where residual non-equilibrium 
effects of AGN ionisation enhance metal columns. 

The general behaviour that we find can be summarised 
as follows. The enhanced photo-ionisation from the AGN 
rapidly raises the ionisation states of metals, thus allowing 
high-ionisation species to exist at higher densities than they 
would in equilibrium with the normal EGB. After the AGN 
turns off, the long recombination times of metals, along with 
the multi-ionisation state nature of metals and their sigiufi- 
cant ionisation fractions (i.e. fx- > 0.1 for multiple i), mean 
that metal columns often only return to equilibrium on a 
timescale that is longer than the AGN-on phase. Thus, a 
whole new range of ionisation models that have not been 
considered are allowed. 



3.1 AGN fossil zones at 2 = 0.9 

In our fiducial example, we consider the effects of an en- 
hanced ionisation in the proximity zone around an AGN that 
is on for 20 Myr at z — 0.9. This case is directly applicable 
to observations by COS at a redshift where a large array 
of both low- and high-ionisation species can be observed. 
This example is motivated by the extremely st rong Nevill 
(A^NcViii ~ lO^^cm"^) absorber observed by I Tri pp et al.l 
(2011) at 2 = 0.927 along the PG1206+459 sight line with 
well-aligned lower ionisation species. A foreground galaxy, 
labelled 1779 by that work, found at the same redshift at 
an impact parameter of 68 kpc, shows signs of being a post- 
starburst galaxy with a B-band absolute magnitude of —22 
and even AGN signatures from [Nev] emission that likely 
cannot be explained by star formation alone. 

Figure [5] shows the time evolution of gas enriched 
to solar metallicity with tih = 10~'*cm~"^, which corre- 
sponds to an overdensity of 77 at 2 — 0.9 (assuming 
Qb = 0.046 and a hydrogen mass fraction of 0.75). This 
choice represents a typical density of metal-enriched CGM 
gas. The gas is assumed to lie at its equilibrium temper- 
ature, T = 10"-^^ K, when we mimic the instantaneous 
turn-on of the AGN by adding our template AGN spec- 
trum to the HMOl EGB at that redshift. We assume an 
AGN flux of /ll = 10"^**-^ergs"^ cm"^ Rz'^ , which has the 
Hi ionisation equivalence of a uniform EGB with JLL.cqu ~ 
■|^g-20.0gj.g g-i pj^-2 jjj^gg].a^g(j over A-k steradi- 

ans; this corresponds to 44 x the EGB value and a total 
ionisation parameter log[?7]= 0.3, where U is the ratio of hy- 
drogen ionising photon density over tih. Using our template 
spectrum and integrating over the soft X-ray band (0.5-2.0 
keV), we calculate that this is the fiux at 100 kpc from a 
io. 5-2.0 ~ lO''*'^^ erg s^^ AGN, which appears to be a rea- 
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sonable assumption given the observed impact parameter of 
68 kpc for galaxy 1779 . 

In the top panel we show the column densities of sev- 
eral observationally accessible ions by assuming an absorber 
length of 52 kpc, which is the predicted size of an Hi ab- 
sorber at this de nsity and temperature according to Equa- 
tion 3 of ISchav"3 ([2001,). All column densities scale linearly 
with the assumed absorber size. As the AGN turns on at 
t — 0, the Hi column density and ionisation fraction drop 
nearly instantaneously by a factor of 44 as Tion.Hi is only 
1100 yr for this fieldQ The ionisation timescales for the metal 
species, shown in the middle panel of Figure [2] during the 
AGN-on phase (thick dotted lines) are longer, 0.1 — 5 Myr, 
but still less than the assumed AGN lifetime of 20 Myr. Dur- 
ing the AGN-on phase, the species Ovi, Nevill, and Mgx 
first rise to a peak and then fall as they are first ionised from 
lower ions to these states and then to higher ions, owing to 
the large increase in the radiation field. On the other hand. 
Civ, which initially comprises 27% of carbon, decreases dra- 
matically as carbon is ionised all the way to C vil in under 
10 Myr. 

The column densities in the AGN-on phase are ap- 
plicable to proximate quasar absorption line (QAL) sys- 
tems observed near the background source quasar, or to 
intervening QAL systems where the local ionising field is 
significantly enhanced by foreground AGN near the sight 
line to the background source quasar. Absorption systems 
with log[?7]= 0.3 will have Nevill column densities that 
are 80 — 300 x stronger than for the HMOl field where 
log[[/]= —1.3 at nn ~ 10~*cm~'^, while A^'hi will be re- 
duced by a factor of 50, and A^civ will be reduced even 
more. A'^ovi can either increase or decrease depending on 
the AGN lifetime. Mgx is the most enhanced photo-ionised 
species shown, increasing more than 10'' x in strength to 



The temperature increases by 29% to T 



10" 



K 



due to photo-heating by the strongly increased ionising field 
(bottom panel). The temperature increase results primar- 
ily from the ejection of highly energetic electrons by X-ray 
photons leading to high-ionisation metal species (e.g. C v, 
Ovil, and O vill). For this solar-enriched case, metal photo- 
heating dominates over that from Hell, which has a lower 
ionisation potential than metal ions. The photo-heating of 
H I does not contribute to the temperature increase, because 
in ionisation equilibrium the increase in the radiation field 
is exactly compensated by the corresponding decrease in the 
neutral fraction, and the latter is too small to result in sig- 
nificant non-equilibrium heating. 

The AGN proximity zone fossil phase begins after the 
AGN turns ofl^ and the ionisation field switches back to the 
uniform HMOl EGB. The solid thin lines in the middle panel 
indicate the recombination timescales of the various species. 
Although Tree, hi: ^ 10^ yr, A^hi reaches equilibrium in a 
much smaller timescale, ~ 10'*'^ yr, because the equilibrium 
neutral fraction it has to return to is only /hi = 10"'*'^. 
In contrast, fossil ionised bubbles around dead quasars can 
persist for ~ rrcc,Hii during reionisation since the equilib- 



^ Note that the ionisation timescale is shorter than the light 
travel times from the AGN to the absorber, which we do not 
consider. 



rium state is neutral (|Furlanetto et al1l2008l ): however, we 
are considering only cases where an EGB ionises the IGM 
to a high level. A^hi remains slightly smaller than before 
t = 0, because photo-heating increased the temperature of 
the gas, thereby reducing the recombination rate. The cool- 
ing timescale of ~ 1 Gyr means that the gas will take of order 
this time to approach the original equilibrium temperature 
and Hi fraction. 

Even though the Helium-like metal ion recombination 
timescales (solid lines in the middle panel of Figure [Sjl are 
~ 20 — 180 X shorter than for Hll, their recombination lags 
are 100 — 900 x longer than for Hi, because the equilibrium 
ion fractions are a significant fraction of unity. A case in 
point is the timescale of ~ 50 Myr that it takes C iv columns 
to approach their previous equilibrium values. A series of 
recombination timescales of species between Cvil and Cv 
determine the total time it takes C iv to return to equilib- 
rium, since carbon is primarily in C vil when the AGN turns 
off. However, the longest is Trcc,cv ~ 40 Myr, shown as the 
green solid line in the middle panel. 

O VI is significantly enhanced during ~ 70 Myr after 
the AGN turns off, even though its column density is de- 
creased during the AGN-on phase. Ovi achieves at least 
twice as large columns between 33-74 Myr, as oxygen re- 
combines from a majority of O vill at t = 20 Myr to O iv 
at t = 95 Myr. Trccovii = 20 Myr is the bottleneck slowing 
the recombination sequence (orange solid line in the middle 
panel) . 

The most dramatic residual observable change is for 
Nevill, for which columns increase 100 fold for the first 15 
Myr of the fossil phase and remain at least 10 x stronger for 
50 Myr. Despite the main recombination bottleneck being 
T"rcc,Noix ~ 10 Myr, or about half as long as rioc,ovi, Nevill 
takes several recombination times to approach its very low 
equilibrium value. Thus A^NeViii ~ 10^'*cm~^, which has 
been observed and attributed to coUisionally ionised gas 
at T ~ 10*^ K ISavage et al.1 l2005l : [Narayanan et al.l |2009| . 
l201lD . could have an alternative colder origin if nearby 
galaxies currently or recentl y harboured an AGN. In passing, 
we note the recent work of Muzahid et al.l l|2013l ). who find 
extremely strong high ionisation lines (e.g. Nevill, Mgx) in 
associated QSO absorbers, which are possibly photo-ionised 
by the AGN; however, we note that those absorbers are 
probably located much nearer to the supermassive black hole 
than the extended proximity zones that are the focus of the 
present study. 

3.1.1 Parameter variations 

We have considered the very specific example of an AGN 
proximity zone fossil applicable to enriched gas at CGM 
densities (i.e. overdensities > 50). While our choices appear 
reasonable, we wish to broaden the range of physical param- 
eters and situations by varying the metallicity, density, AGN 
strength, and AGN lifetime. We will relate these variations 
to situations that may occur in the real Universe and will 
argue that fossil zones could significantly change the inter- 
pretation of metal QAL observations. We will continue to 
assume that the metals are cool (T ^ 10'* — 10*'^ K), re- 
siding at the equilibrium temperature where photo-heating 
balances radiative cooling. 

Metallicity: Choosing solar metallicity for an absorption 
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Figure 2. The evolution of a njj = 10~^cm~'^, 52 kpc, solar abundance absorber irradiated by the z = 0.9 HMOl EGB when an AGN 
with intensity JLL.oqu = 10~^''ergs~^ cm~^ Hz~^ sr~^ (^0.5— 2.0 = 10**'^ergs~^ at 100 kpc distance) turns on at t = for 20 Myr. After 
20 Myr the AGN turns off but the long recombination times of Helium-like metal ions (C V, O VII, Ne IX, &; Mg XI; between 5-50 Myr; 
middle panel, solid thin lines) cause these observable Lithium-like ions to take longer than the AGN-on time to reach their equilibrium 
levels. lonisation timescales of the Lithium-like ions (middle panel, thick dotted lines shown only when the AGN is on) are shorter than 
20 Myr when the AGN is on. Photo-heating, primarily of metals, results in a slightly increased temperature of the gas that cools on a 
timescale of 1 Gyr (middle panel, thick black dashed lines). Unlike metals and despite a 1 Gyr recombination timescale, hydrogen rapidly 
recombines to equilibrium when the AGN turns off, owing to its very small equilibrium neutral fraction. The columns scale linearly with 
the assumed absorber length and metallicity. 
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system is not unprecedented as detailed analysis of ob- 
served system s often yield so l ution s with solar or supe r -solar 
metallicities (iPeroux et all l2006l: JProchaska et alj l2006l : 



ISchave et al]|2007l : lTripp et al.ll201l| - |Tumlinson et alll201lh 

However, many of these systems are attributed to lower 
density gas, which may on average be enriched to lower 
levels if metallicity increases with density as observa- 
tions and simu lation s indi c ate (e 
Pavel I2OO6I , I2OO9I : 



Schave et al 



Wiersma et al 



200a 



201c 



Qppenheimer 

2OIII ). Lowering the metallicity does not alter the recom- 



bination timescales or behaviour significantly. The equilib- 
rium temperature in the HMOl ionisation field increases 
from r = 10*-^^ K for Z = Z© to lO''''^ K for 
0.1 Zq, owing to reduced metal-line cooling, and this re- 
sults in only a 0.2 dex reduction of A'hi- Hence, metal- 
line strengths as a function of time can, to first order, 
be scaled with Z. However, if intergalactic metals are 
poorly mixed on s mall scales as some observations and 
simulations indicate (Simcoe et al. '200d';' Scha ve et al. II2OO7I: 



2011]), 



lOppenhc imcr & Dave 2009; Tcppcr-Garcia et al.l 
then high metallicities may even be appropriate at low den- 
sities. 

Density: Changing the density alters the non-equilibrium 
behaviour because the recombination times scale as nu~^- 
The ionisation parameter also changes when assuming a con- 
stant radiation field, meaning that in equilibrium atomic 
species are photo-ionised to higher ion states at lower den- 
sities. We explore here four densities: nu = 10"^", 10~*'^, 
10~* °, and 10~"^'^cm~^ corresponding to z — 0.90 overden- 
sities p/{p) = 7.7, 24, 77, and 243 respective ly. We assume 
an absorber length set by the ISchavd l|200ll ) relation {I = 
0.52nH~'^'^\/T/10* kpc), which equals the Jeans scale in a 
gas of uniform density and temperature. We set T — 10* K 
in the previous equation for simplicity, even though tempera- 
tures are usually slightly higher. This size scale, which ranges 
from 164 to 29 kpc for our densities, has been show n to re- 
produce t he observed sizes o f the Lyg abso rbers (e.g |Schav3 
I2OOII : [Prochaska ct al. 2004: lLehner et al.L 2007') and agrees 
very well with predictions from hydrodynamic simulations 
JPave et al.ll20ld : iTepper- Garcia et al.|[20l3 : iRahmati etlH 
|2012| ). 

Figure [3l top panel, shows that the A^hi of the ab- 
sorbers initially range from 10^^'^ — >■ lO^^'^cm"^ with in- 
creasing density, which correspond to observed Lya ab- 
sorbers t hat typicall y exhibit de t ectable metal-l i ne ab- 
sorption JPan forth fc Shulll l2008l: iThom fc Cheiil l2008bl : 
iTripp et al.1l2008l : ITilton ct al. 2012). As the AGN radiation 
field turns on and off, A'^ei reduces and increases by a factor 
50 on timescales that are too short to appear in the figure. 
A slight deviation occurs for higher densities where photo- 
heating raises the temperature during the AGN more, owing 
to the greater number of bound electrons available to ionise 
from metals. Although lower densities have a smaller photo- 
heating temperature increases (bottom panel), the longer 
cooling times of several Gyr mean that a 20 Myr AGN phase 
has longer lasting effects on the temperature of the IGM. 

The third panel in Figure [3] shows A^ovi for the four 
densities. The overdensity range we explore covers nearly the 
same range that simulations pr edict observed O vi absorbers 
to arise from at z < 0.5 (e.g. [O ppenheimcr & Dave 2009; 
Tepper-Garcia et al.ll2011 ; Cen fc Chisari 2011; Smith et al.l 
I2OI2I '). WhenThe AGN turns on. 
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Figure 3. The density dependence of the non-equilibrium 
proximity effect at z = 0.9 using an AGN with Jll.bqu = 
-|^Q-20.0gj.gg-l (,j-|^-2 jj^-l ^^urning on at t = for 20 Myr. 
The assumed absorber lengths decrease from 164 to 29 kpc as 
density increases from njj = 10"^" to 10~^'^cm~^. A^hi, A^civ, 
A^OVli A^NcVllli ^iid temperature evolution are shown from top 
to bottom. Ionisation timescales remain unaltered between the 
four cases, but recombination times scale as njj"^ meaning that 
the fossil effect will last longer at lower density. Metallicity is 
assumed to be solar in all cases, but metal columns scale with 
metallicity with the minor exception that photo-heating will be 
lower at lower metallicity and temperatures will not vary as much. 
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the ionisation timescales are independent of riH, but the 
new ionisation equihbrium does depend on nn. For nn ^ 
10~''cm~^ oxygen is ionised from O vi to higher states re- 
ducing these columns, while for nn ^ 10~^'^cm~^ lower 
ionisation oxygen is ionised to O vi increasing the columns. 
When the AGN turns off, lower densities take longer to re- 
turn to ionisation equilibrium because recombination times 
scale inversely with nu, Thus, while the lowest density will 
have deficient O vi relative to the equilibrium case for > 100 
Myr, the highest overdensity will have enhanced O vi for a 
shorter timescale, similar to the AGN lifetime of 20 Myr in 
this case. 

In contrast to O vi, A'^civ is reduced by the AGN for 
all densities while A'^ncViii shows an increase for nn S5 
10"* '^ cm"^ (2nd & 4th panels of Figure [3] respectively). The 
C IV ionisation fraction peaks at nu = 10~'^'^cm~^ when 
photo-ionised by the z = 0.9 HMOl field at T ~ 10* K. 
Turning on the AGN shifts this peak to nu = 10~^'^cm~^. 
Hence, C iv will be reduced by the AGN, and the lowest 
densities will take the longest to recombine from higher 
ions in the fossil zone. Conversely, the ionisation fraction 
of Neviii peaks at nu = lO'^^'^cm"^ at T ~ 10* K, 
which shifts to 10"^-^cm~^ with the AGN. Hence, gas with 



10^ 



2x108 



10 



"^■^cm ^ can have TVngViii columns as high as 



lO^'^cm"^ during the AGN phase, 12,000x larger than for 
the HMOl field in equilibrium. The recombination times of a 
few Myr means that this strong Ne vill will be shorter lived 
in the fossil zone. 

AGN radiation field increase: In Figure |4l we con- 
sider field increases, JLL.equ of factors 4.4 and 440 to sim- 
ulate the AGN proximity effect at 3.16 and 0.32 x the 
fiducial radius (i.e. 316 and 32 kpc) respectively for the 
same strength AGN. Or alternatively, an AGN 0.1 and 
10 X as strong at 100 kpc. We again use our fiducial den- 
sity nn = 10~* °cm~^. The ionisation timescales decrease 
in proportion to the field increase, and stronger radiation 
fields cause metals to become more ionised in a shorter 
amount of time. However, post-AGN metal-line strengths 
show qualitatively similar evolution, because the recombi- 
nation timescales are the same. There is a greater delay for 
stronger fields, refiecting the additional recombination time 
from even higher ionisation states (e.g. Trcc.cvi, tvcc.oviii). 



Even for Jll,( 



10 erg s cm Hz sr , which cor- 



responds only to a factor 4.4 x increase relative to the EGB, 
significant enhancements are seen, especially in Ne Vlll. The 
field increase is insufficient to ionise to Ne ix, but post-AGN 
Nevill columns are enhanced 100 x and decline slowly ac- 
cording to Tree, Ne VIII. Considering that the 10 x smaller field 
increase applies to 10^''^ ~ 32x more volume, the weaker 
field case could be most relevant to intervening QAL, which 
provide a volume-weighted sampling of the IGM. 
AGN lifetime: Thus far we have assumed a single tagn = 
20 Myr lifetime, but AGN lifetimes are predicted to be as 
short as 10® yr. Figure[S]shows what happens if an AGN is on 
for shorter periods at multiple times. We choose to show an 
AGN that is on 5 out of every 50 Myr yielding a duty cycle, 
d = 10%, and an AGN that is on 1 out of every 25 Myr, 
yielding d = 4%. For metals, the effect of a shorter lifetime 
is similar as that of a weaker AGN: they are not as highly 
ionised by the AGN, but the post-AGN recombination times 
are the same. A shorter tagn results in metals not reaching 
as highly ionised states. Post-AGN, the metals are less likely 
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Figure 4. The dependence of the non-equilibrium proximity ef- 
fect on AGN strength at z = 0.9 with nu = 10~''cm~3 and I = 52 
kpc for an AGN turning on at t = for 20 Myr. The AGN flux en- 
hancement relative to the HMOl EGB ranges from a factor of 4.4 
(long dashed lines) to 440 (short dashed lines). The stronger the 
AGN flux, the higher the ionisation state that the metals reach, 
and the more delayed the recombination to equilibrium becomes 
owing to the need to recombine through more ionisation states 
with significant recombination times. Although the metallicity is 
assumed to be solar, the metal columns scale with metallicity, 
except for the minor effect of photo-heating being less efficient at 
lower metallicity. 
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to have to recombine through these upper states reducing 
the number of recombination times to return to equiUbrium. 
Nevertheless, the most significant bottleneck is the longest 
metal Tree from the Helium-like state (Cv, Ovii, Neix) to 
the Lithium-like ions we plot. Hence metal-line evolution is 
often similar in these three cases although shifted in time 
depending on the strength and of course the frequency of 
AGN-on phases. 

Even short AGN lifetimes can efficiently enhance metal 
columns if the duty cycle is similar to or smaller than the 
recombination times. A case in point is Ne vill for tagn = 1 
Myr in the 4th panel, whose column density is enhanced by 
60 X averaged over 0-200 Myr, despite the AGN being on 
only 4% of the time. For tagn = 1 Myr, O VI is enhanced 
by a factor of three and remains relatively constant, owing 
to a coincidence: Tion.ovi ~ 1 Myr is similar to tagn, and 
Troc.ovii ~ 20 Myr is similar to the duty cycle timescale of 25 
Myr. The implications for small duty cycles with resonances 
that amplify the fossil effect will be explored further in i]3.3l 
Finally, the bottom panel shows how duty cycles can pump 
temperatures up, owing to the combination of short photo- 
heating timescales and long cooling timescales if the gas has 
solar metallicity. 

Extra-galactic background and AGN EUV power 
laws: We have chosen an A GN spectrum with an EUV 
power law oi j^-^-^"^ based on iTelfer et~ai] (|2002l '). The lo- 
cal AGN ionises diffuse gas that is additionally ionised by 
the HMOl spectrum, which itself assumes a QSO EUV in - 
p ut of If we used ei ther the lHaardt fc Madaul (| 19961 ') 
or lHaardt fc Madaul (|2012l ') EGBs that use hard QSO input 
spectra of and u~^'^'^ respectively, then the increase 
in the radiation field at the ionisation potentials of high 
ions like Ovi, Nevill, and Mgx would be more muted. The 
general behaviour explored above would still be applicable, 
but the ionisation equilibrium with the EGB would be more 
highly ionised. 

In the real Universe, there exists an extremely large 
ran ge of observe d AGN EUV power laws. Jclfcr ct all (|2002l ') 
and IShull et al.l (|2012l ) together find a range oi v ^ to 1^° 
encompassing most of their measured AGN, and even an 
outlier with i^"'^®, which is the well-known quasar HE2347- 
4342 at z = 2.9 that has b een used as a probe of Hen 
reionisation (|Shull et al.ll2O10l ). Hence, the fossil effect may 
be even greater for harder AGN, which are more capable 
of ionising high ions, whi le softer AGN will have much less 
of an effect. Interestinglv. IScott et al.] (|2004l 'l provide a sam- 
ple showing an anti-correlation between EUV hardness and 
AGN luminosity. Their composite s pectrum h as a harder 
EUV slope, u-°-^'^, than lTelfer et all (2002) and lShull etHI 
(|2012l ), which suggests greater fossil effects. In fact, they 
find that their lowest luminosity AGN have EUV slopes ex- 
ceeding 0.0 in the low-redshift Universe, which makes our 
exploration of weak AGN in tj3.3l possibly even more rele- 
vant. 



3.1.2 Direct modelling of the z = 0.927 PG1206+459 
TVevill absorber 

We show more detailed results for our fiducial model in Fig- 
ure[Hl where we plot the evolution of the column densit ies of 
some metal- line species observed bv I Tripp et al.l (|201ll ) and 
H I against Ne vill for the example shown in Figure [S] We 
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Figure 5. Different AGN lifetime and duty cycle variations 
of the non-equilibrium proximity effect at at 2 = 0.9 with 
njj = 10~^cm~^ and I = 52 kpc for an AGN flux of JLL.oqu = 
10~'^''erg cm"'^ Hz~^ sr"'^. We compare our fiducial case of a 
single AGN-on phase of 20 Myr with an AGN-on phase of 5 Myr 
every 50 Myr and 1 Myr every 25 Myr. The first AGN on phase 
always starts at t = 0. Significant non-equilibrium behaviour oc- 
curs in all cases significantly altering the time-averaged column 
densities during the AGN-off phase. Photo-heating in bursts can 
boost average gas temperatures as well. Although the metallicity 
is assumed to be solar, the metal columns scale with metallicity, 
except for the minor effect of photo-heating being less efficient at 
lower metallicity. 
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plot the coloured diamonds every 10^ yr starting at t = 
when the AGN turns on (large magenta diamond). The up- 
per left panel shows the evolution of Nv and Nevill. The 
ionisation timescales imply that these ions take several Myr 
to approach the AGN-on equilibrium. At 20 Myr the AGN 
turns off (indicated as a white line in the colour bar in the 
lower left panel), and the evolution back to the normal HMOl 
EGB equilibrium is much slower, taking over 1 00 My r. 

The open squares show the iTripp et al.l (|2011 ) com- 
ponent absorbers, indicating lower limits for saturated ab- 
sorbers as triangles in the case of O iv, and upper limits as 
upside-down triangles for non-detections of Mgx. We com- 
bine their -f65 and -|-108kms^^ components into one moti- 
vated by the blended Ne vill detection. While photo-ionised 
Nevill is far too weak without an AGN, the AGN ionisa- 
tion achieves the observed Nevill columns within 0.5 Myr. 
A possible scenario for this absorber is an AGN that has 
been off for ~ 30 Myr, as most Neviii, Nv, Sv, Mgx, and 
O IV components are then within the allowed ranges. How- 
ever, Mgll is underestimated by up to a factor lO'*, and Hi 
is too low even after the AGN has turned off. This illustrates 
the need for a multi-phase model with a component for the 
lower ionisation species that is shielded from the AGN or 
which arises in much denser gas. 

Another possible scenario in Figure |S] is the AGN just 
turning on for a few x 10^ yr producing reasonable columns 
for higher ionisation species (Nv, Sv, Neviii, Mgx), but 
requiring a second phase for lower ionisation species (Hi, 
Mgll, Oiv). After aU, galaxy ITTg does show active AGN 
signatures. This scenario is, however, a priori less likely since 
10^ < 10^ yr. 

Figure |S] shows merely one possible time-dependent his- 
tory for an AGN ionising cool, metal- enriched gas. A weaker 
AGN strength provides nearly as good a fit. A higher den- 
sity, riH = 10~^'^cm~^, model where the AGN turned off 10 
Myr before fits the columns just as well, but produces 100 x 
stronger Mgll, which is however still ^ 100 x too low. Our 
point is that models with time-dependent AGN that have 
either turned off in the past, recently turned on, or fluctu- 
ate are worth considering, especially for the extreme case of 
the PG1206-I-459 absorber with evidence of a post-starburst 
AGN within 100 kpc of the sight line. The difRcuhy is that 
a wide range of models are possible, most of which are un- 
constrainable without knowing significantly more about the 
galaxy-absorber pair. 

We emphasise that a proximity zone fossil can enhance 
metal lines over a significant path length if the AGN was 
strong in the past. We showed that a wide range of AGN 
intensities, as would be expected along a large path length 
through a proximity zone, will lead to similar recombina- 
tion time histories and hence similar enhanced columns for 
high-ionisation species at a given density. If a proximity 
zone is many Mpc in size, then enhanced metal columns 
should be observed along a similar path length. A fossil zone 
could be traced by unusually strong high-ionisation species 
ove r a large vel o city r ange as in the case of the Nevill in 
the iTripp et al.l (|201ll ) absorber, which extends over 1500 
kms"'^. This corresponds to a Hubble fiow of over 10 proper 
Mpc implying an extended fossil zone and a very strong 
AGN in the past. 

In summary, we propose a new time-dependent model 
for the z = 0.927 PG1206-f459 absorber that relies on en- 



hanced photo-ionisation of T ~ 10* K gas from either an 
AGN that turned off tens of Myr ago or an AGN that 
recently turned on. While our model can reproduce the 
columns of many of the higher ionisation lines, including 
Nevill, it requires a multi-phase structure with absorption 
by lower ionisation states, especially Mgll arising in much 
denser gas remaining shielded from the AGN. This scenario 
pr ovides an alt e rnativ e to the outflow shell model proposed 
by iTripp et all (|201ll ) of a T > lO'^ K hot gas component 
and a velocity-coincident cold gas component, which may 
be expected for a propagating shock with post-shock radia- 
tively cooling metal-enriched gas. Our AGN model requires 
a much smaller mass for high-ionisation components, as lit- 
tle as 3 X lO^'M© of gas, assuming a path through the cen- 
tre of a spherical r ~ 14.5 kpc, solar-enriched cloud with 
riH ~ lO-^ '^cm-^ 



3.2 AGN and proximity zone fossils at the peak 
of QSO activity at z ^ 2.5 

Another application of our code is the transverse proximity 
effect (TPE) where a foreground QSO ionises gas along the 
sight line towa rd a background QSO . Recent work has been 
)erformed bv iHennawi et al.l (|2006|). [ Worseck fc Wisotzkil 
20061 ): I Worseck et all l|200i1 ), and iGoncalves et al.l (|2008l ) 
who explored the TPE at z ~ 2.5, near the peak of QSO 
activity. The TPE can be elusive when traced by Hi alone 
owing to complicating effects of the overdense regions cre- 
ating more Hi absorption, which goes in the opposite direc- 
tion of decreas ed Hi owing to the Q SO ionisation (see also 
ISchirber et aLli2004 ). For example, iHennawi fc Prochaskal 
l|2007h argue that QSO emission is anisotropic based on 
the enhancement of Lyman-limit systems in transverse sight 
lines indicating that the QSO radiation does not effect off- 
sight line directions as often, even after accounting for the 
fact that QSOs may l ive in overdense regions. 

IGoncalves et al.l l|2008h analysed metal-line systems in 
a quasar triplet at z 2.5 and found that the TPE was 
apparent when C iv and O vi aligned with H I are included 
in the analysis. They in fact argued that unusual Novi/Nui 
and A^civ /Nm ratios of components in transverse proximity 
zones indicate that the QSO radiation field reaches off-sight 
line directions, and that the two foreground QSOs, which 
are extremely bright and rare, more likely radiate isotropi- 
cally. Their work shows that the non-linear ionisation prop- 
erties of metal lines, when used i n conjunction with Hi, can 
help constrain the T PE (see also lWorseck fc Wisotzki|[2006l : 
IWorsecket aLlbOOTh . 



3.2.1 The QSO-on phase 

We turn on our AGN template spectrum for 20 Myr and 
find the average change in column density owing to the en- 
hanced field. The left tables in Figure [7] show, from top 
to bottom, the average enhancement in Hi, Civ, Nv, and 
Ovi of the time-averaged column densities during the 20 
Myr QSO-on phase for a grid of densities (log[nH/cm~'^]) 
and fluxes (log[JLL,equ/(erg s"'^ cm~^ Hz~^ sr~^)]). The sec- 
ond column in each table shows equilibrium column densi- 
ties (log[iV/cm~^]) using the HMO l z = 2.5 EGB, absorber 
lengths according to lSchav3 l|200ll l. and solar metallicity. As 
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Figure 6. Column densities of observed species in tiie PG1206+459 z = 0.927 absorber by iTripp et al] l|201ll l (black squares are 
observed detections, triangles lower limits, upside down triangles upper limits) compared to the evolution of column densities predicted 
for the case shown in Figure [2] njj = lO^^cm""^, T = lO* '^ K, Z = 52 kpc, solar metallicity with an AGN with intensity JLL,oqu = 
lO^^^ergs"^ cm~^ Hz~^ sr~^ cm"'^ turning on at t = and off at t = 20 Myr after which the radiation field returns to the ever present 
2 = 0.9 HMOl EGB. Nv, Sv, Mgx, Mgll, Oiv, and Hi are plotted against Nevill in different panels from upper left to bottom right. 
Time t = is shown as the large magenta diamond, and small diamonds are plotted every 10^ yr, colour-coded by age with the AGN-on 
phase being coloured magenta. The white line in the colour bar indicates when the AGN turns off and is obvious in the inflection point 
in the N V and Mg X plots. This model can approximately reproduce observations of Ne VIII, N V, S V, Mg X, and O IV at ~ 50 Myr (30 
Myr after the AGN turns off), but severely underestimates Mgll and also tll. 



with all cases in this paper, we assume the gas initially rests 
at the temperatures where photo-heating balances cooling 
for the HMOl field, which fall from lO"*-^** ^ 10* °^ K as 
density increases from nn — 10~^ 10"'^ cm"'' at z — 2.5. 
The grid cells to the right of the HMOl column indicate the 
column density change (51og[Af] (in dex) relative to the HMOl 
equilibrium, colour-coded to be redder if the QSO decreases 
the column and bluer if it is increased. 

The behaviour of Hi is relatively easy to understand. 
As the AGN turns on, A^hi equilibrates to a column density 
that is inversely proportional to the intensity of the new 
radiation field. We define the AGN enhancement factor as 

JeGB + J'AGN,cqu /„^ 
SAGN = z . (6) 

Jegb 

With the z = 2.5 HMOl EGB having log[JLLl=-21.3 com- 
bined with the weakest AGN field log[JLL,cqu]=— 21.5, Hi is 
reduced by log[5rAGN]= 0.2 dex. The six AGN field enhance- 
ments we list range from pagn=1.6-200, which decrease A'hi 
in the same proportion modulo the effect of photo-heating 
raising the temperature of the gas and thus further reducing 
Hi. Photo-heating by the strongest AGN reduces Hi by an 
additional factor of two in nn — 10~^cm~'' gas, owing to 
the ejection of energetic metal electrons at X-ray energies. 



The behaviour of the metal lines is not as simple as 
that of Hi. As we saw in the previous subsection, AGN al- 
most always ionise C iv to higher states and thus reduce 
their columns, while enhancing O vi at higher tih and re- 
ducing it at lower nn, depending on the field strength. We 
again assumed solar metallicity, but the listed 51og[A'^] values 
are applicable to other metallicitics modulo the effect of in- 
creased photo-heating at higher metallicity. The same trends 
are observed bv lGoncalves et al.l (|2008l ): increased Ovi and 
decreased C iv. Their qagn ranged from 10 — 200 and they 
derived densities of tih = lO^'''^ — 10~^'^cm~^, which would 
correspond to the bottom three rows and right-most four 
columns in the left tables of Figure [T] One difference is that 
they use a single power law of in their UV EGB, while 

we are using the HMOl field enhanced with a AGN 
spectrum plus an X-ray component. This allows our field 
to ionise O vi to O vil and above, which explains why we 
find more reduction in A^'ovi for the stronger fields: Ovi is 
ionised to higher states much more easily. 

An additional point about the TPE during the AGN- 
on phase is that O vi increases significantly more than H I 
decreases at higher densities and lower field strengths (see 
blue cells in bottom left panel of Figure [7]). We may thus ex- 
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Figure 7. Time-averaged column densities of Hi (top panels), Civ (2nd panels), Nv (3rd panels) and Ovi (bottom panels) for a z = 2.5 
AGN-on phase lasting 20 Myr (left panels) and the fossil effect of the 20 Myr immediately after the AGN turns off (right panels). The 2nd 
column lists log[iV/cm~^] for the normal HMOl EGB at 2 = 0.25, assuming ionisation equilibrium and solar metallicity, for 5 different 
densities (listed as log[(cm~^)] in 1st column). The remaining columns indicate the effect of the AGN intensities (listed in the 2nd row in 
log[JLL,equ] units, where log[ JLL,equ]= —21.3 at 2 = 2.5 for the HMOl EGB) shown by the increase or decrease in log[A'^] relative to the 
HMOl log[Af]. Cell colour indicates whether the AGN increased (bluer) or decreased (redder) the column density. Hi is heavily affected 
during the AGN-on phase, but returns to nearly the original column afterwards. C IV, N V, and O VI are also strongly affected during the 
AGN-on phase, and show time-averaged columns over the next 20 Myr that remain out of equilibrium. Solar metallicity is assumed, but 
the metal columns scale with metallicity, except for the minor effect of AGN photo-heating raising the temperature slightly. 



pect significantly enhanced O vi with weaker than usual H I 
toward the edge of a proximity zone, which may approach 
~ 10 proper Mpc in size for the brightest QSOs. In fact, the 
large enhan c emen t in proximate O vi absorbers observed by 
iTripp et al.1 l|2008l ) within 1500 kms"^ of the brightest QSOs 
at lower redshifts (or ~ 17 Mpc at z = 0.5), combined with 
their lower observed Nui/Novi ratios agrees with this pre- 
diction. We put forth a model that these are not absorbers 



directly associated with the QSOs, but instead intervening 
IGM absorbers on the outskirts of the proximity zone with a 
Qagn of a few an perhaps even less than two. This of course 
assumes that the AGN has a harder EUV spectrum than 
the EGB. 

A similar point can be made that N V could be enhanced 
in AGN proximity zones. Only cells with nn = 10~^ ''cm~'^ 
and weaker AGN fields show noticeably enhanced N v when 
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the AGN is on, but these ceUs may be well-represented in 
QSO proximity zones. Consider that while N-^y for the 
equilibrium HMOl field peaks at jih ~ 10~"^'^cm~ (assum- 
ing Zq as in Figure[7]), the highest observed N v column may 
correspond to riH = 10~'^" given that metallicity likely cor- 
relates with density. Factoring in that most of the proximate 
volume and path length corresponds to the weaker fields, the 
strongest observable Nv columns could correspond to the 
outskirts of proximity zones. Since nitrogen is less abundant 
than other metals and N v is much rarer than C iv and O vi, 
this AGN enhancement could boost normally unobservable 
or hard-to-detect Nv columns by ~ 5 — 10 x into obvious 
detections. 



3.2.2 QSO proximity zone fossils 

We now cross over to the fossil phase by considering the 
time-averaged columns of Hi and metals during the first 
20 Myr after the QSO turns off. This timescale is chosen 
to demonstrate that significant fossil effects can persist on 
a similar timescale as tagn. We show the deviations in 
columns on the right side of Figure [T] Listed is the time- 
averaged column deviation over this time span (note that 
the actual deviation will likely be larger or smaller at a spe- 
cific fossil zone time). 

The Civ column is almost always reduced, but not to 
the same extent as during the QSO-on phase. Nonetheless, if 
a significant fraction of intervening C iv absorbers reside in 
fossil zones, then models th at try to fi t these assumi ng a uni- 
form EGB like HMOl fe.g.lSchaye et al. 2003: Simcoe etal] 



Vagn, exposed to at least qagn — 1 is 



l2004l : IOppenheimer fc Davell2006l : lTescari et aLH201li ') would 
underestimate the amount of carbon. Hence, the true global 
ionisation correction of C iv would be larger, and the in- 
ferred IGM metallicity would be too low. The opposite would 
be true for Ovi at higher densities, but it would have the 
same trend as C iv if it mostly arises from lower density gas. 
Meanwhile, Hi would remain nearly unaffected except for 
the slightly reduced columns owing to photo-heating of met- 
als. The resulting absorption line systems will have ratios 
(e.g. Ovi/Hl and Civ/Hl) that would differ significantly 
from those predicted by equilibrium models with a normal 
EGB or in a proximity zone. 



3.2.3 How affected are metal lines by fluctuating AGN? 

We now present an estimate of the fraction of the volume of 
the Universe that may be affected by transverse AGN prox- 
imity and fossil zones at high-z. We estimate the luminosity 
function ($(L^, #/Mpc~Vlog(Li. l) of AGN at 912 A us- 
ing the analytic fitting formulae of iHopkins et al.) (|2007l ) at 
z = 2.a For a given isolated AGN, the proper radius within 
which the AGN fiux is at least qagn — 1 times is given by 



nosity at v in ergs s 



1) J^), where is the AGN lumi- 
Hz~^. Hence, the fractional volume. 



Vagn 



0(L, 



47r 



3 \{4nyigA 



1)J. 



3/2 



We used their script at 
|http://www. cfa.harvard.edu/~phopkins/Site/qIf. html[ 



{l + zfdL, 
(7) 

when integrating over the luminosity function {u is in the 
rest frame). The (1 -I- z)^ term is a result of $(ii^) being in 
comoving Mpc and our proximity zone units being proper 
units. Using J„ = 10^^^'^erg s~^ cm"^ Hz~^ sr~^ at 2 = 2.5, 
the volume ionised by AGN to at least the level of the HMOl 
background, qagn ~ 2, is 0.16%, and the volume ionised to 
<?AGN = 10 is 0.005%. This calculation does not depend on 
how anisotropically AGN emit, because the average sub- 
tended solid angle over which AGN of a given luminosity 
emit is compensated by the number of all emitting AGN for 
that luminosity. The Vagn fractions are much lower at lower 
redshift, where AGN are less frequent and the Universe has 
expanded. 

While the fraction of volume with qagn ^ 2 is small 
compared to the volume filling factors of metals in simu- 
lations calibrated to fit obser ved Civ statistics, which is 
of order 10% at z ~ 2 — 3 (jOppenheimer fc Pavel [2OO6I : 
iBooth et aLll2012l ). there are reasons to argue that AGN 
proximity and fossil zones could affect a significant frac- 
tion of observed metal lines. First, strong C iv systems 
{Nciv lO'^^cm"'') are associated with Ly man-break galax- 
ies (LBGs ) at small radii (80-85 proper kpc lAdelberger et al.l 
2005; Stci del et al.ll20l"c)l ). and the co-moving line density of 
such strong absorb ers is consistent with being within 50 kpc 
of ^ 0.5L* LBGs (jCooksev et al.ll2012l ). This implies that 
many of the strongest metal-line absorbers occupy the small 
fraction of volume with enhanced ionisation. 

More importantly, the multiplicative effects of AGN 
with short duty cycles can significantly expand the volume 
affected by fossil zones. Our lower redshift example in Fig- 
ure [5] shows that a d = 4% duty cycle leads to 3 x greater 
O VI and 5 x reduced C iv columns time-averaged once the 
AGN first turns on at t = 0. This duty cycle is relevant 
since t he observed fractio n of LBGs with active AGN is 
3% l|Steidel et al.ll2002l ), which could refiect their duty 
cycle. Taking 1/d for d — 4% means that the fossil effect 
applies to a 25 x greater volume than the proximity zone 
effect, and hence that 4% of the total volume is affected by 
5AGN ^ 2. This volume fraction is similar to that thought 
to be requir ed to account for the observed incidence of weak 
metal lines (jOppenheimer fc Davell2006l : iBooth et al.ll2012l ) 
suggesting that non-equilibrium AGN fossil zones may be 
critically important. Finally, metal lines can be enhanced if 
<?AGN < 2, as is shown in Figure[7]for log[JLL,Gqu]=— 21.5, or 
5AGN = 1.6. This leads to a doubling of the volume affected 
by both proximity and fossil zones, and Civ is reduced by 
2 — 4x in almost every case. Again, this arises because the 
AGN EUV spectrum is harder than that of the EGB. 

We conclude that a large fraction of metal lines could 
arise from fossil zones at high-z, but precisely what fraction 
depends on many factors. If the AGN-on phase is long-lived 
and singular for most AGN then the fossil effect would be 
limited. However, if AGN are fluctuating with short duty cy- 
cles in typical star-forming galaxy hosts and have hard EUV 
slopes, then the fossil effect could be a critical consideration 
for typical metal lines. We may have to re-evaluate the phys- 
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Figure 8. The time evolution of O VI columns (top panel) and gas 
temperature (bottom panel) for a npj = 10~^cm~^, I = 52 kpc 
solar metallicity absorber at z = 0.25, wiiere an AGN of the listed 
intensity turns on for 1 Myr every 10 Myr starting at t = 0. The 
periodic AGN activity significantly enhances O VI columns and 
doubles the gas temperature over 200 Myr in the most extreme 
case via photo- heating of metals. This case could be applicable 
for metals around normal looking galaxies if AGN duty cycles are 
short and have frequencies comparable to the O VI recombination 
timescale. 

ical interpretation of Civ and Ovi absorbers at high- 2:, and 
the amount of underlying metals traced in diffuse gas. 



3.3 Low-redshift Seyferts ionising their 
circumgalactic media 

Our last case study are low-redshift metal absorbers asso- 
ciated with weaker AGN, such as Seyferts. We consider an 
admittedly idealised situation of a 10% duty cycle where the 
AGN turns on for 1 Myr out of every 10 Myr at z = 0.25. We 
plot the time history of Ovi in Figure [8] for three different 
field enhancements, to show how the continuous pumping 
resulting from a short duty cycle can even significantly en- 
hance the ionisation around weak AGNs. We again assume 
riH ~ lO^^cm"'^ solar-enriched gas at the equilibrium tem- 
perature of T = 10"-^ K. 

This situation may be applicable to Markarian galaxies 
in the local Universe. These are among the biggest, brightest 
spirals, and have an unusually blue central core that usually 
requires non-thermal emission, which often indi cates AGN 
activ ity. About 10% of Markarians are Seyferts jWeedmanl 
1 19771 ). which motivates our choice of a 10% duty cycle. Of 
course, not all Markarians are Seyferts with a short duty cy- 
cle, but we consider the possibility that ~10% of Markarians 
were more active AGNs in the recent past (i.e. in the last 
~ 10 Myr). The three cases plotted in Figure [HI JLL.cqu = 
10~^^ °, 10-21°, and 10-2° °ergs-i cm'^ Hz^^ sr'^ corre- 
spond to soft X-ray luminosities at 100 proper kpc of 



1/0.5-2.0 = \ 10^'^ \ and 10'**-'ergs-i AGN, respec- 
tively, using our AGN template spectrum. These luminosi- 
tie s are in the range of Seyfert X-ray luminosities observed 
bv lRush et al.l l| 19961 ). The bottom panel of Figure [8] shows 
that the temperature doubles to 10*'* K in the strongest 
case. 

We tabulate the time-averaged enhancements of C iv, 
Ovi, and Nevill over 200 Myr, or the first 20 cycles, 
for a grid of densities (log[nH/cm-'^]) and AGN fluxes 
(log[JLL,oqu/(ergs-i cm-2 Ifz-i sr^i)]) in Figure (9] Similar 
to Figure [T] the second column in each table shows col- 
umn densities {\og[N /cm~'^]) using the H MOl EGB (her e 

at 2 = 0.25), absorber lengths according to ISchav3 (|200lf ). 
T = 103-85-4-40 ^ |3gfQ^g ^-jjg ^^j.jjg ojj^ g^^^ g^jg^j. metal- 
licity. The grid cells to the right indicate the column density 
change (51og[A'^] relative to the HMOl equilibrium, colour- 
coded to be redder if the column density goes down due to 
the AGN and bluer if it goes up. The wh densities corre- 
spond to overdensities ranging from 27 to 2700. While we 
show solar metallicity, the change is similar for lower metal- 
licities. 

The main point of this diagram is to show that even for 
a 10% AGN duty cycle, the time-averaged columns are sig- 
nificantly changed, mostly in a positive direction for these 
higher metal ions. While there has never been a survey de- 
signed to probe the halos of galaxies showing signs of on- 
going or recent AGN activity in the low-redshift Universe, 
it may be worth explorin g in light of the COS observations 
of iTumlinson et all ((2011) revealing very strong O vi around 
normal spiral, star-forming galaxies at 2 ~ 0.2. Markarian 
galaxies, for example, may show even stronger Ovi if this 
gas originates from a photo-ionised component at T < 10^ 
K. 



4 SUMMARY 

The assumption of a uniform extra-galactic background 
breaks down in an AGN proximity zone, which is over- 
ionised relative to the rest of the IGM. Hi columns are 
reduced and metals are ionised to higher states owing to 
increased ionisation. We demonstrate here that once the 
AGN turns off, metals residing in "fossil" AGN proximity 
zones may be over-ionised for timescales exceeding the typ- 
ical AGN -on phase. Using the non-equi librium code intro- 
duced in lOppenheimer fc Schavd l|2013h . we follow metal- 
enriched gases in itially ionised by the ex tra-galactic ionisa- 
tion background l|Haardt fc Madciull200ll ) under the effect of 
increased ionisation from a local AGN. After the AGN turns 
off, metal ions in fossil proximity zones often require many 
Myr to return to the equilibrium values, owing to significant 
recombination times of metal ion species. We demonstrate 
how this non-equilibrium fossil effect alters observational 
metal-line diagnostics used to constrain the physical state 
of intergalactic and circum-galactic gas. In particular, we 
show that commonly observed Lithium-like ions (i.e. C iv, 
Nv, Ovi, Neviii) can be heavily affected at a range of red- 
shifts from the present-day to z^ 3. This behaviour is in 
contrast to Hi, which re-equilibrates comparatively rapidly 
once an AGN turns off, owing to its equilibrium state having 
very low neutral fractions in the IGM. 

We consider metal-enriched gas with densities nu ~ 
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Figure 9. Time-averaged column densities of Civ (top panel), 
Ovi (middle panel), and Nevill (lower panel) assuming that a 
z = 0.25 AGN turns on for 1 Myr every 10 Myr. The 2nd col- 
umn lists log[Ar (cm-2)] for the normal HMOl field at z = 0.25, 
assuming solar metallicity, for 5 different nn densities (listed as 
log[(cm~^)] in the 1st column). The rest of the columns indicate 
the time averaged effect of the AGN intensities (listed in the 2nd 
row in log[JLL,cqu] where log[ JLL,oqu]= —22.2 at z = 0.25 for 
the IfMOl EGB) shown by the relative increase or decrease in 
log[A''] relative to the HMOl log[A''] . Cell colour indicates whether 
the AGN increases (blue) or decreases (red) the column density. 
Despite the AGN being on only 10% of the time, metal columns 
can be affected, mostly in a positive direction. Solar metallicity is 
assumed, but the metal columns scale with metallicity, except for 
the minor effect of AGN photo-heating raising the temperature 
slightly. 

10"^ - 10"^cm"^ aX z = 2.5,0.9 and 0.25, which corre- 
sponds to diffuse IGM and CGM gas, that is initially in 
thermal equilibrium (T ~ 10* — 10'''^ K). Our parameter 
choices are based on hydrodynamic simulations showing a 
significant fraction of diffuse metals in thermal equilibrium, 
where cooling balances photo-heating, owing to short cooling 
times from higher temperatures (T ~ 10^ — 10^ K). The gas 
begins in ionisation equilibrium with the HMOl EGB, and 
becomes irradiated by an AGN template spectrum turned 
on for 1 to 20 Myr corresponding to an ionisation increase 
of between 1.6 — 500x the EGB radiation field strength at 
the Lyman limit. 

Our fiducial case corresponds to a local AGN turning 



on at 2; = 0.9 where a wide range of metal lines are observ- 
able with current facilitie s (e.g. COS on HST) and inspired 
by the lTripp etHI l|201lD observations of very strong Ne vill 
absorption in a system associated with a galaxy exhibiting 
post-starburst and possible AGN signatures. We show that 
for reasonable AGN energies, Nevill columns are increased 
by a factor of ~ 100 during a 20 Myr AGN-on phase and 
also the 15 Myr in the following fossil phase, and remain 
> 10 X stronger for 50 Myr after. We show for a variety 
of gas densities, AGN strengths, and AGN lifetimes that 
Nevill columns of 10^'* — lO^^cm"'' are possible for signifi- 
cant timescales during the fossil phase fo r solar abundances . 
We directly model the components of the lTripp et al.l (|201ll ) 
absorber, finding that a model with T < 10^ K gas in a fos- 
sil zone of 10-30 Myr can reproduce many of the observed 
columns and upper limits using as little as 3 X lO^'M© of 
solar enriched gas per component. Solutions are also possi- 
ble if the AGN is currently on. However, a separate phase is 
necessary for all solutions to explain the Mg ll and strong H I 
columns. Our application to this absorption system demon- 
strates that a whole new class of non-equilibrium solutions 
are available for metal-enriched diffuse gas in the presence 
of evolving ionisation fields. 

We also model AGN ionisation at z ~ 2.5, including 
the transverse proximity effect for a foreground QSO zone 
in a QSO pair, where observed Civ/Hl and Ovi/Hl line 
ratios indicate enhanc ed ionisation over the normal EGB 
(jGoncalves et al.]|2008l ). We first show that active proxim- 
ity zones can exhibit these unique ratios, and then consider 
how the non-equilibrium fossil effect can alter C iv and O vi 
columns for timescales longer than the AGN lifetime. Unique 
line ratios are possible where H I has rapidly re-equilibrated 
but the metal lines are still recombining. We argue that it is 
possible that fossil zones could affect many more absorbers 
in the high-z Universe than AGN proximity zones if AGNs 
fluctuate with short duty cycles, have hard extreme UV 
slopes, and are associated with typical star-forming galax- 
ies, which prolifically enrich the IGM. A large fraction of 
intervening metal absorbers in typical QSO spectra could 
be affected by the fossil effect, meaning that ionisation cor- 
rections assuming a uniform EGB would be wrong. C iv 
columns are almost always reduced in fossils meaning that 
the true diffuse gas metallicity becomes higher than when 
one assumes equilibrium with the EGB field. 

Finally, we consider lower redshift AGN ionising their 
surrounding CGM with a 10% duty cycle where the AGN is 
on for 1 out of every 10 Myr. Despite this short duty cycle, 
high ions (Ovi and Nevill) are significantly enhanced of- 
ten by more than 1 dex, and temperatures are increased 
by a factor of two for stronger fields. The key to these 
increases of metal columns and temperature are that the 
photo-ionisation/photo-heating timescales are short or com- 
parable to the AGN lifetime, and the recombination/cooling 
timescale is long or comparable to the duty cycle interval. 
The result is intriguing in that surprisingly little energy in- 
put is necessary to greatly enhance metal columns in the 
halo of fluctuating typical AGN. 

These non-equilibrium fossil and proximity zone effects 
may complicate the analysis of intergalactic metal lines that 
assume photo-ionisation equilibrium with a uniform EGB. 
The fossil effect is unique in that H I re-equilibrates rapidly 
while metal lines remain over-ionised leading to different 
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ionisation corrections to obtain metallicities. This could be 
particularly important during the quasar era, z ^ 2 — 5, for 
which we find that the fraction of the volume occupied by 
fossil zones may be similar to the fraction that is thought to 
be enriched with heavy elements. 
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